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Towards the feasibility study of the electroweak baryogenesis in realistic technicolor scenario, we 
investigate the phase structure of 2+iVf -flavor QCD, where the mass of two flavors is fixed to a small 
value and the others are heavy. For the baryogenesis, an appearance of a first order phase transition 
at finite temperature is a necessary condition. Using a set of configurations of two-flavor lattice QCD 
and applying the reweighting method, the effective potential defined by the probability distribution 
function of the plaquette is calculated in the presence of additional many heavy flavors. Through the 
shape of the effective potential, we determine the critical mass of heavy flavors separating the first 
order and crossover regions and find it to become larger with Nf. We moreover study the critical 
line at finite density and the first order region is found to become wider as increasing the chemical 
potential. Possible applications to real 2 + 1-flavor QCD are discussed. 
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Introduction. Understanding the phase structures 
of QCD at the finite temperature and density is an in- 
teresting subject. Numerical simulations based on lattice 
gauge theory have played an essential role there. Espe- 
cially, the quark mass dependence of the nature of the 
phase transition in 2 + 1-flavor QCD is one of the impor- 
tant topics in this research field. The phase transition is 
first order when all of the up, down and strange quark 
masses are sufficiently small or large, and is crossover 
in the intermediate region between them. Starting from 
the massless three-flavor QCD and increasing the strange 
quark mass only, chiral transition is expected to change 
from first to second order at a tricritical point . 

The techniques developed in numerical simulations of 
lattice gauge theory works effectively to study a class 
of new physics models beyond the standard model (SM) 
as well. Technicolor (TC) model [l| built of many fla- 
vor QCD (or equivalently vector-like SU(3) gauge the- 
ory with many fermions transforming as the fundamental 
representation) is one such model. In this class of models, 
the Higgs sector is replaced by a new strongly interacting 
gauge theory and its spontaneous chiral symmetry break- 
ing (SxSB) causes electroweak (EW) symmetry breaking. 
Among various attractive features, we focus on the pos- 
sibility of the EW baryogenesis in this scenario Q , which 
requires TC gauge theories to go through a strong first 
order chiral phase transition. As mentioned for the QCD 
case above, the nature of the phase transition depends 
on the number of flavors and masses 0] . In realistic TC 
models, two flavors of them are exactly massless and the 
mass of other Nf flavors must be larger than an appro- 
priate lower bound otherwise SxSB produces too many 
(light pseudo) Nambu-Goldstone (NG) bosons. Three 



of them are absorbed into the longitudinal mode of the 
weak gauge bosons, but any other NG bosons have not 
been observed yet. On the other hand, the first order 
transition at small mass terminates at the critical mass 
like 2 + 1-flavor QCD. Thus, if one requires the first order 
EW phase transition in TC model, it brings in the up- 
per bound on the mass of Nf flavors. This motivates us 
to study the thermal nature of 2+Af-flavor QCD. As we 
discuss below, the critical mass increases with Nf . Hence, 
the boundary of the first order region can be investigated 
more easily for large Nf. 

Another purpose of this study is to understand the 
real QCD with 2+1 flavors. At the physical masses of 
up, down and strange quarks, the chiral transition is 
crossover at zero density, and is expected to become first 
order at a critical density. Finding this critical density is 
one of the most interesting topics in the lattice QCD sim- 
ulations and in the heavy-ion collision experiments. To 
this end, the simulation research of the critical surface 
in the masses and chemical potential parameter space is 
important 0, @j. However, recent lattice QCD studies 
suggest that the critical region at zero density is accessi- 
ble only when the quark masses are very small and thus 
its determination is difficult Q. The study of 2 + Nf- 
flavor QCD is a good testing ground for investigating 
iVf-indepcndent universal properties, such as the critical 
scaling near the tricritical point. This will provide im- 
portant information for 2 + 1-flavor QCD. 

Method. To study the phase transition, we calculate 
the effective potential defined by the probability distri- 
bution function of the plaquette. The distribution func- 
tion has two peaks at a first order transition, since two 
phases coexist with the same probability. The nature of 
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the transition can be thus identified through the shape of 
the potential @, 0] ■ We define the plaquette distribution 
function for 2+A/f-flavor QCD with the quark masses m,f 
and chemical potential [if (/ = !,■■• , Nf ) by 



w{P;P,m f ,n f ) = VUV^V^S(P-P) 



~S Q S a 



N { +2 

VU 5{P-P) e 6/3Ar - teP I (detM(m/,/i/)),(l) 



where S g and S q are the gauge and quark action, respec- 
tively, and M the quark matrix. A s i te = A s 3 x N t is the 
number of sites. (3 = is the lattice bare parameter. 
P is the generalized plaquette operator, and this method 
is applicable to the case of improved actions replacing 
P to P = —S g /(6N S itc[3). Normalizing by the partition 
function, Z — J w(P)dP, cq. ([T]) gives the histogram for 
P. The effective potential is then given by 



V eS {P;P,m f ,[j, f ) = -]nw(P;l3,mf,fif). 



(2) 



We consider QCD with two degenerate light quarks 
of the mass mi and Nf heavy quarks with a mass m^, 
and the chemical potential fi is imposed only on the two 
light quarks, for simplification. Denoting the potential 
of two-flavor QCD at zero density by Vb(P;/3), that of 
2+iVf-flavor QCD is written as 

V eS (P;l3,m h ,v) = Vo(P-,i3o)-lnR(P;p,m il ,[i;j3o), (3) 
with 



lnR(P;p,m h ,ii;po) = 6(/3 - (3 )N site P 



6{P - P) 



det M{m h /i) VYdet M (m h , 0) x N 
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det M (mi, 0) J \ det M(oo, 0) 



(4) 
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where (•••} l g denotes the ensemble average over two- 
flavor configurations generated at /?, m\ and the vanish- 
ing [i. Since the m\ dependence is not discussed in the 
following, mi is omitted from the arguments. /3q is the 
simulation point, which may differ from j3 in this method. 
By performing simulations at various /?o, one can obtain 
the potential in a wide range of P. 

Restricting the calculation to the heavy quark region, 
the second determinant for Nf flavors in eq. ([5]) is ap- 
proximated to the leading order as 



In 



detM(« h ) 



= 288iV B i to KjP + l2N 6 s {2K h ) Nt VL +• • -(5) 



det M(0) 

for the standard Wilson quark action and 
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(6) 



' S(P- P)exp[6hN^Q}' 

R(P; h, 0) = 

S(P' - P] 



for the four-flavor standard staggered quark, in cq. (O 
is the hopping parameter being proportional to 1/rah, 
and VL is the real part of the Polyakov loop. For improved 
gauge actions such as S g = 6iV sito /3[c (plaquette) + 
ci (rectangle)] , additional c\ x 0(k 4 ) terms must be con- 
tained in eqs. (|5|) and ©, where c\ is the improvement 
coefficient and cq = 1 — 8ci. However, since the improve- 
ment term does not affect the physics, we will cancel 
these terms by a shift of the coefficient C\ . 

At a first order transition point, V e ff shows a double- 
well shape as a function of P, and equivalcntly the cur- 
vature of the potential d 2 V e ff/d 2 P takes a negative value 
in a region of P. To observe this behavior, /3 must be 
adjusted to be the first order transition point. However, 
from cqs. ^ and ([5]), d 2 V c ff/dP 2 is independent of /?. 
The fine tuning is not necessary in this case Q. More- 
over, d 2 V ff/dP 2 over the wide range of P can be easily 
obtained by combining data obtained at different (3. Wc 
therefore focus on the curvature of the effective potential 
to identify the nature of the phase transition. 

Denoting h = 2Nf(2n\ l ) Nt for the Wilson quark ac- 
tion, or h = Nf/(A x (2m-f 1 ) Nt ) for the staggered quark 
action, we obtain In R(P; 0, Kh, 0; 0o) = \nR(P;h,0) + 
(plaquette term) + 0(n^ t+2 ) for fj, = with 



(7) 



Notice that R(P; h 7 0) does not depend on (3 . The pla- 
quette term does not contribute to d 2 V c ff/dP 2 and can 
be absorbed by shifting fi -> (3* = f3 + 48N { k^ for Wilson 
quarks. As seen above, the choice of the quark action 
is not important. In the following, we discuss the mass 
dependence of R through the parameter h. 

Numerical results. We use the two-flavor QCD 
configurations generated with p4-improved staggered 
quark and Symanzik-improved gauge actions in Ref. 0], 
thus P = —Sg/(6N sitc (3). The lattice size A^ito is 
16 3 x 4. The data are obtained at sixteen values of j3 
from (3 = 3.52 to 4.00 keeping the bare quark mass to 
ma = 0.1. The number of trajectories is 10,000 - 40,000, 
depending on j3. The corresponding temperature normal- 
ized by the pseudo-critical temperature is in the range 
of T/T c = 0.76 to 1.98, and the pseudo-critical point is 
about j3 — 3.65, where the ratio of pseudo-scalar and 
vector meson masses is mps/mv ~ 0.7. All configura- 
tions are used for the analysis at zero density, while the 
finite density analysis is performed every 10 trajectories. 
Further details on the simulation parameters are given in 
Ref. @. The same data set is used to study the phase 
structure of two-flavor QCD at finite density in Ref. Q. 

We first calculate the potential in two-flavor QCD at 
zero density, Vo(P;(3), the first term in eq. ([3|). Be- 
cause the finite temperature transition is crossover for 
two-flavor QCD at a finite quark mass, the distribution 
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FIG. 1: Top: In R(P; h,0) as functions of the plaquette. Bot- 
tom: The curvature of hxR(P;h,0) for h = 0.01 - 0.07. The 
circle and square symbols are d 2 Vo/dP 2 (P). 



function is always Gaussian type. We thus evaluate the 
curvature of Vq using an identity for the Gaussian distri- 
bution, d 2 Vo/dP 2 = 6N sitc /xp, where xp is the plaque- 
tte susceptibility, X p = 6N aite ((P - (P)) 2 ). The slope of 
Vq in the heavy quark limit can be also measured using 
an equation derived from eqs. ((3]) and ([5]). When one per- 
forms a simulation at /3o, the slope is zero at the minimum 
of Vo(P; /3o), and the minimum is realized at P sa (P)f3 - 
Hence, we obtain dV Q /dP((P) ft) , p) = -6(fi - l3 Q )N sitc 
[Tfjj . The result of d 2 V /dP 2 is plotted in the bottom 
panel of Fig. [TJ The circle symbols with dashed lines 
are calculated by xp- The square symbols are computed 
by the numerical deferential of dVo/dP obtained at the 
minimum of V$. dVo/dP are the squares in Fig. [2] These 
results obtained by two different methods are consistent. 

Zero density. In the calculation of R(P]h,0), we 
use the delta function approximated by S(x) ~ 
l/(AV7r)exp[-(a;/A) 2 ], where A = 0.0025 is adopted 
consulting the resolution and the statistical error. Be- 
cause R(P; h, 0) is independent of (3, we mix all data 
obtained at different j3 as is done in Ref. Q ■ The results 
for In R(P; h, 0) are shown by solid lines in the top panel 
of Fig. [T] for h = 0.01 - 0.07. A rapid increase is observed 
around P ~ 0.82. It is also important to note that the 
gradient becomes larger as h increases. 

The second derivative d 2 \nR/dP 2 is calculated by fit- 
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FIG. 2: The mass dependence of the slope of V e ff(-P; /3, h, 0) 
normalized at (/3, h) — (3.65, 0) for h = 0.0 - 0.1. The squares 
are dV /dP. 



ting In R to a quadratic function of P with a range of 
P ± 0.015 and repeating with various P. The results are 
plotted in Fig. [T] (bottom), where d 2 Vo/dP 2 is also shown 
as the circles or the squares with dashed lines. This fig- 
ure shows that d 2 (In R)/dP 2 becomes larger with h, and 
the maximum around P = 0.81 exceeds d 2 Vo/dP 2 for 
h > 0.06. This indicates that the curvature of the effec- 
tive potential, d 2 V cS /dP 2 = d 2 V /dP 2 - d 2 (ln R)/dP 2 , 
vanishes at h ~ 0.06 and a region of P where the curva- 
ture is negative appears for large h. We estimated the 
critical value h c at which the minimum of d 2 V c g/dP 2 
vanishes and obtained h c = 0.0614(69). 

To see the appearance of the first order transition in a 
different way, we plot dV e g/dP at finite h for j3* = 3.65 
in Fig. [21 The shape of the dV s/dP is independent of /3 
because d 2 V e s/dP 2 is /3-independent. dV e ff/dP is mono- 
tonically increasing when h is small, indicating that the 
transition is crossover. However, the shape of dV a s/dP 
turns into an S-shapcd function at h ~ 0.06, correspond- 
ing to the double-well potential. 

We defined the parameter h = 2N{ x (2nh) Nt for the 
Wilson quark. Then, the critical Khc corresponding h c 
decreases as Khc = [h c /(2N{)] 1 ^ Nt /2 with N{, and the 
truncation error from the higher order terms in be- 
comes smaller as Nf increases. The application range 
of the hopping parameter expansion was examined in 
quenched QCD simulations with N t = 4, by explicitly 
measuring the size of the next-to-leading order (NLO) 
terms of the expansion [ll[. They found that the NLO 
contribution becomes comparable to that in the leading 
order at Kh ~ 0.18. Hence, this method may be applica- 
ble up to around ~ 0.1. For instance, in the case of 
Nf = 10 with Nt = 4, K hc is 0.118. 

Nonzero density. Finally, we turn on a chemical po- 
tential yu for two light quarks and discuss the fi depen- 
dence of the critical mass. We keep the chemical potential 
zero for the other quarks. Because the strange chemical 
potential is known to be small in the heavy-ion collisions, 
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this situation for 2 + 1-flavor corresponds to the experi- 
ments. As discussed above, we can investigate the criti- 
cal region easily for large Nf. R(P; h, fi) is then given by 
the expectation value of (det M{m\, fx)/ det M(m\, 0)) 2 x 
(detM(m h ,0)/detAf(oo,0)) Wf with fixed P. The quark 
determinant is computed using the Taylor expansion of 
ln[detM(mi,/x)/detM(mi,0)] in terms of fx/T up to 
0[(fi/T) 6 ) and the Gaussian approximation is applied to 
avoid the sign problem as explained in Ref. Q ■ This ap- 
proximation is valid for small fx. The truncation error has 
been estimated comparing the results up to 0(fx 4 ) and 
0{fi 6 ) for fi/T < 2.5 and is found to be small_0. The left 
panel of Fig. [3] shows the curvatures of In R(P; h, fx) at 
fx/T = 1 and Vq. The maximum value of the curvature of 
In R(P; h, fi) is larger than that at fx = 0. This means the 
critical h is smaller at finite \i. Figure [3] (right) shows the 
critical value of h as a function of fx. In the region above 
this line, the effective potential has the negative curva- 
ture region, indicating the transition is of first order. It 
is clear that the first order region becomes wider as fx 
increases. If the same behavior is observed in 2 + I-flavor 
QCD, this gives the strong evidence for the existence of 
the critical point at finite density in the real world. 

Although this analysis is valid only for large Nf, it 
gives a frame of reference for the study of critical mass 
at finite fi. Notice that lnR(P; h, fx) is given by the 
sum of In R(P; 0, fx) and In R(P; h, 0) approximately and 
that the behavior of ln^(P; h,0) in Fig. [T]is very sim- 
ilar^ to that of lnR(P;0,fx) in Figs. 5 and 7 in Ref. Q. 
lni?(P; 0, fx) is estimated from the quark number suscep- 
tibility at small fi and In R(P; h, 0) is obtained from the 
Polyakov loop at small k^. Both the quark number sus- 
ceptibility and the Polyakov loop rapidly increase at the 
same value of P near the transition point, which enhances 
the curvature of \nR. Therefore, the critical h decreases 
with fx or equivalently the critical fx decreases with h. 
The same argument is possible for 2 + 1-flavor. 

Conclusion and outlook. We studied the phase 
structure of 2+Af-flavor QCD to explore the realizabil- 
ity of the EW baryogenesis in technicolor scenario and to 
understand properties of the finite density QCD. With 
fixing the mass of two light quarks, we determined the 
critical mass of the other Nf quarks separating the first 
order and crossover regions. The critical mass becomes 
larger with Nf. Furthermore, the chemical potential de- 
pendence of the critical mass is investigated for large Nf, 
and the critical mass is found to increase with fi. 

To quantify the strength of the first order phase tran- 
sition is the next important step for the estimation of 
the baryon number asymmetry in TC scenario. An- 
other interesting application of our method is to study 
universal scaling behavior near the tricritical point. If 
the chiral phase transition in the two flavor massless 
limit is of second order, the boundary of the first or- 
der transition region mf(mh) is expected to behave as 
mf ~ |TO^ rl - — m/J 5 / 2 in the vicinity of the tricritical point , 
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FIG. 3: Left: The curvature of In R(P; h, fi) as functions of 
the plaquette at fi/T = 1.0. Right: The critical line in the 
(h, fi) plane. In the region above this line, the transition is 
of first order. The square at fx = is computed using all 
configurations, and the others are measured with every 10. 



(mi,mh,/i) = (0, m^ rl , 0), from the mean field analysis. 
This power behavior is universal for any Nf. The density 
dependence is important as well, which is expected to be 
m \ ~ I A* 1 5 [S- Starting from large Nf, the systematic 
study of properties of QCD phase transition is possible. 
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